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ABSTRACT 

Recently, Shen et al. (2005) studied the contributions of curvature effect of 
fireballs to the spectral lag and showed that the observed lags could be ac- 
counted for by the effect. Here we check their results by performing a more 



2 ' precise calculation with both the formulas presented in Shen et al. (2005) and 

^ ' Qin et al. (2004). Several other aspects which were not considered in Shen 

et al. (2005) are investigated. We find that in the case of ultra-relativistic 
, motions, both formulas are identical as long as the whole fireball surface is 

■ concerned. In our analysis, the previous conclusion that the detected spectral 

lags could be accounted for by the curvature effect is confirmed, while the 
conclusion that the lag has no dependence on the radius of fireballs is not 
true. We find that introducing extreme physical parameters is not the only 
outlet to explain those observed large lags. Even for the larger lags (~ 5s) (see 
Norris et al. 2005), a wider local pulse {Ate,FWHM — W^s) could account for 
it. Some conclusions not presented in Shen et al. (2005) or those modified in 
our analysis are listed below: a) lag cx F"*^ with e > 2; b) lag is proportional 
to the local pulse width and the FWHM of the observed light curves; c) a 
large lag requires a large ao and a small Pq as well as a large E^p] d) when 
the rest frame spectrum varies with time, the lag would become larger; e) 
lag decreases with the increasing of Re, f) lag oc E within the certain energy 
range for a given Lorentz factor; g) lag is proportional to the opening angle 
of uniform jets when 9j < O.GF^^. 
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1 INTRODUCTION 



It is common that there is a spectral lag between different energy channels in gamma-ray 
bursts (see, e.g. Cheng et al. 1995; Norris et al. 1996, 2000; Wu & Fenimore 2000; Chen 
et al. 2005). Although There are several attempts to explain the origin of the spectral lag 
(See, e.g., Salmonson 2000; loka & Nakamura 2001; Schaefer 2004; Kocevski et al. 2003a), 
the problem remain unresolved. 

Recently, Shen et al. (2005) (hereafter Paper I) tentatively studied the contribution of 
curvature effect of fireball to the lag, and the resulting lags are very closed to the observed 
one. In fact, in their paper introducing extreme physical parameters such as F < 50 or 
a > —0.5 is not the only outlet to explain those observed large lags. A reasonable intrinsic 
pulse width t'^, i.e., hundreds to thousands of seconds, which could produce the lag of 0.1 — 1 
s based on their model, could also account for it. 

Qin (2002) had derived in a much detail the fiux intensity based on the model of highly 
symmetric expanding fireballs, where the Doppler effect of the expanding fireball surface is 
the key factor to be concerned, and no terms are omitted in his derivation, therefore the 
formula is applicable to the cases of relativistic, sub-relativistic, and non-relativistic motions 
and to those of spherical fireballs or uniform jets. With this formula, Qin (2003, 2004) derived 
the observed photon count rate of a fireball (See Eq. (1) in the paper). 

A detailed comparison between Shen et al. (2005) and Qin et al. (2004) model (hereafter 
Paper II) gives rise to the conclusion that the two models are identical when the following 
conditions are satisfied. Firstly, as shown in equation (8) of Paper II, the two quantities, 
local time tg and observed time t, are related hy t — {l — (3cos9){tg — tc)+tc+D/c—RcCos9/c, 
where tc is the initial local time, Rc is the radius of the fireball measured at tg = tc, and D is 
the distance of the fireball to the observer. When taking tc = and T = t + (Rc — D)/c, one 
would find that the above relation will become the one determined by equation (3) in Paper I 
(where T is referred to as observed time defined in Paper I). Secondly, in equation (3) of Paper 
II, the integral lower and upper limit are determined by ^min = cos'^ min{cos^inin, {te,max — 
t + D/c)/{{te,m3^-tc)P + Rc/c)} and ^^ax = cos'^ max{cos9^^, {t0,^in-t + D/c)/{{te,nim- 
tc)P + Rc/c)}, respectively, where ^min and ^max are confined by the concerned area of the 
fireball surface, and te.min and te.max are confined by the local time range of emission. For 
such a local pulse, te^min = and t^^max oo, when taking 6'min = and 9ma.x ~ 7i'/2 (i.e., 
the photons reaching the observer come from the whole surface of the fireball), one would 



obtain 6mm = and ^max = cos^^{l — T/t) (where r = Rc/c), which are consistent with 
the integral lower and upper limit in equation (4) of Paper I, respectively. In a word, the 
mainly difference between the two models is that, Shen et al. neglected the constraint of 
the local pulse and the concerned area of fireball surface on the integral lower and upper 
limit of the observed flux (see equation (4) in their paper) , which plays a important role in 
producing the observed light curves. Strictly, when dealing with uniform jets with a small 
opening angle such as 6* < l/F, the formula presented in Shen et al. (2005) would fail to be 
applied. 

We in this paper check their results by performing a more precise calculation with the 
two models, at the same time, several other aspects which were not considered in Shen et al. 
(2005) are investigated, such as lag's dependence on energy and opening angle (in Section 
2). The previous conclusion that the observed spectral lags could be accounted for by the 
curvature effect is confirmed and some new results are obtained. For example, lag decreases 
with the increase of radius of fireball other than weakly dependent on it. Lag oc F"*^, e > 2 
in the case of fixing the local pulse width, which is identical with that, lag ex F~^ in the case 
of fixing the intrinsic pulse width, obtained by Shen et al. (2005). Accordingly an interesting 
question arises: would the two approaches always lead to identical results? which will be 
discussed in Section 3. And then we give our discussions and conclusions in the last section. 



2 THEORETICAL ANALYSIS ON SPECTRAL LAGS 

Under the assumption that a fireball expands isotropically with a constant Lorentz factor 
F > 1 and the radiation is independent of direction, the expected count rate of the fireball 
measured within frequency interval [vi, 1/2] can be calculated with equation (21) in paper II, 
which can be rewritten as follows 

2nRl P'— J(r,)(l + Pre)\l - r + r,)rfr, /- ^-^^dv 

/icL'2F3(l-/3)2(l + ^r)2 ■ ^' 

Note that the formula is in terms of the integral of r^, which is a dimensionless relative local 
time defined by tq = c{t0 — tc)/Rc, where te is the emission time in the observer frame, called 
local time, of photons emitted from the concerned differential surface dse of the fireball {9 is 
the angle to the line of sight), tc is the initial local time which could be assigned to any values 
of tff, and Rc is the radius of the fireball measured at = tc. Variable r is a dimensionless 
relative observation time defined by r = [c{t — tc) — D + Rd/Rc, where D is the distance of 



the fireball to the observer, and t is the observation time measured by the distant observer. 
In formula (1), /(r^i) represents the development of the intensity magnitude of radiation in 
the observer frame, called as a local pulse function, and go,i^{j^Q,e) describes the rest frame 
radiation mechanisms. 

Light curves determined by equation (1) are dependent on the integral limits re^min and 
T^e,ma.x, which are determined by the concerned area of the fireball surface, together with the 
emission ranges of the radiated frequency and the local time. But for commonly adopted 
mechanisms, such as the bremsstrahlung, Comptonized and synchrotron radiations, they 
cover the entire frequency band, and do not provide constraints on the integral limits (see 
equations (14) and (15) in Qin 2002). Thus the integral hmits are only determined by 

Te,min = max{re,min, (t - 1 + C0s6lmax)/(l " /3cOs9uiax)} and fe,max = min{T0,inax, (t - 1 + 

COS 6'inin) / (1 —P COS 6'inin) }, whcrc Tectum and rg^max are the lower and upper limit of tq confining 
I{Tg), and ^^inin and 6ma.x are confined by the concerned area of the fireball surface. The 
radiations are observable within the range of (1 — cos^min) + (1 ~ /9 cos 6*111111)761,111111 < t < 

(1 - cos 6'inax) + (1 - /? COS 6'iiiax)Te,max (sCC Paper II). 

Because peak times of different light curves associated with different frequency intervals 
[i/i, 1/2] or different energy bands [Ei, £'2] are different, following Shen et al. (2005), we define 
the spectral lag as the time between the peaks of the light curves in two different channels, 
a lower energy channel [Ei,E2] and a higher energy channel [-£'3,-^4]. In the following, we 
will investigate on what parameters the lags are dependent in terms of the fireball model. 

For the sake of simplicity, we employ in the following a local Gaussian pulse and employ 
the Band function as a rest frame radiation spectrum which was frequently and rather 
successfully employed to fit the spectra of GRBs (see, e.g., Schaefer et al. 1994; Ford et al. 
1995; Preece et al. 1998, 2000) to study the issue. A local Gaussian pulse and the Band 
function are written as 

I{Tg) = /oexp[-(^^^^)2] ^^^^^^ < (2) 

(7 

and 

(i^MU+ao exp[-(2 + ao)^] (^ < ^) 

' ' (frS)"°-*^-P(/5o-ao)(Sf)^^^° (Sj^f^) 

respectively, where Jq, a and To^min are constants, ao and /5o are the low and high energy 
indices in the rest frame, respectively, and i/o,p is the rest frame peak frequency. Typical 



values of the low and high energy indices coming from statistical analysis are ao = —1 and 
Po = -2.25, respectively (see Preece et al. 1998, 2000). 

Due to the constraint to the lower limit of tq, which is Te^rnin > — (see Paper II), 
we assign tq^ — 10(7 + 

^e,min SO that the interval between To^jnin ^iid to^q would be large 
enough to make the rising phase of the local pulse close to that of the Gaussian pulse. 
Prom (2) one can obtain Ato^fwhm = 2\^ln2a, which leads to u = Atq^fwhm /2\/ln2, 
where Atq^fwhm is the FWHM of the Gaussian pulse. From the relation between Tg and tg, 
one gets Ate,FWHM = {Rc/ c)^'Te,FWHM- In the following, we assign {2tiR^^Iq) / [hcD"^) = 1, 

i?c = 3 X 10^^ cm, Te^rnin = 0, 9min = 0, and 9max = 7r/2. 



2.1 Lag's dependence on the Lorentz factor and the local pulse width 

As shown in Qin et al. (2004), the Lorentz factor and the local pulse width are important 
in producing the light curve observed. We wonder if the Lorentz factor and the local pulse 
width play important roses in producing a lag. The spectral lag of GRBs could be found 
either in hight BATSE energy channels or in lower Ginga ones (see, e.g., Gheng et al. 1995; 
Ghen et al. 2005; Norris et al. 1996, 2000; Wu & Penimore 2000). We consider the following 
three different energy channel pairs: BATSE channels 1 and 3, BATSE channels 1 and 4, 
and Ginga channels between 2-10 keV and 50-100 keV (marked as Lagis, Lagi^ and Lago, 
respectively) . 

As suggested by observation, the value of the peak energy Ep of the GRB's uF^, spectrum 
is mainly distributed within 100 ~ 600 keV (See Preece et al. 2000) and peaks at about 250 
keV. Here we assume 250 keV and 50 keV as the typical values of the peak energy of hard 
and soft bursts, respectively. In addition, we assume that typical hard and soft bursts come 
from the same Lorentz factor of the expanding motion of the fireball surface and differ only 
in their rest fame peak energy Eg^p. As derived in Qin (2002), when taking into account the 
Doppler effect of fireballs, the observed peak energy would be related with the peak energy 
of the typical rest frame Band function spectrum (ckq = —1, Po — —2.25) by Ep ~ 1.67P£'o,p. 
In a recent work, the Lorentz factors of a GRB sample were found to be distributed mainly 
within (100, 400) and to peak at about 200 (Qin et al. 2005a). So, when assigning F = 200 
to be the Lorentz factor of these sources, one would find the typical values of the rest frame 
peak energy £^o,p of hard and soft bursts are 0.75 keV and 0.15 keV respectively. 



Presented in Figs. 1 and 2 are the relationships between the lag and the Lorentz factor 
r as well as between the lag and the local pulse width, Atg^pwHM- 



We find from the left panel in Fig. 1 that each of the three kinds of lag has a peak value. 
The Lorentz factor corresponding to this peak value is denoted by Fp (where subscript "p" 
stands for "peak"). The lags increase with F when F < F^; the contrary is observed when 
F > Fp. For Eq^p — 0.75keV, the lags decrease with the Lorentz factor in the law of lag oc F~^ 
with e > 2 for the three kinds of lag. The figure exhibits Lagc < Lagi^ < Lagi^i within the 
range of F > 100 (note that current GRB models suggest that the gamma-ray photons come 
from a relativistically expanding fireball surface with Lorentz factor F > 100, see Piran 
2005). The result and that obtained by Shen et al. (2005) are identical (note that here we 
fix the local pulse width instead of the co-moving pulse width). One can check that, when 
the local pulse width is large enough, for a fixed rest frame peak energy the observed peak 
energy would increase with F following the law of Ep = 1.67F-Eo,p for the adopted Band 
function spectrum (as shown in Qin 2002). According to the law of Ep = 1.67TEo^p, the lag 
will have a peak value when Ep shifting from a higher energy channel to a lower one, which 
could be observed in Fig. 1. Fig. 1 shows that the relationship between the lag and F would 
differ according to different rest frame peak energy Eq^p, which suggests that, in producing 
a lag, not only the Lorentz factor is important but also the rest frame peak energy plays a 
role. 



Fig. 2 shows that the lag increases with the local pulse width following the law of Lag oc 
Ate^FWHM, and it obeys the law of Lagc < Lagi^ < Lagi4, the same as that shown in Fig. 
1. However, as suggested by Fig. 2, the dependence of the lag on the local pulse width is 
much obvious than that on the Lorentz factor F and the rest frame peak energy Eq^p (see 
Fig. 1). In a recent work, Lu et al. (2005a) found that FWHM oc A.to^FWHM based on the 
same model (Qin 2002, Qin et al. 2004) , where FWHM is the full width at half-maximum 
of the GRB's light curve. Thus one could come to lag oc FWHM, which is in agreement 
with the observation that wider pulses were found to produce longer lags (see Norris et al. 
1996, 2001, 2005). 
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Figure 1. -Relationships between Lagis, Lagu, Lagc and F for Eo,p=0.75 keV (the left panel), and that between Lagis and 
r for different Eo,p (the right panel), where the Band function rest frame radiation form with ao = —1 and /3o = —2.25 is 
adopted and we take Rc = 3x 10^^ cm and Atg^pwHM = lO^s. The open circle, the cross and the open square in the left panel 
stand for Lagis, Lagii and Laga, respectively. In the right panel, solid lines from the top to the bottom stand for Bo,p=0.75, 
0.15, 0.075 keV, respectively. 



2.2 Lag's dependence on spectral parameters of the Band function 

We know that different values of the spectral parameters (i.e., ckq) and £^o,p) of the Band 
function are associated with different mechanisms (see Band et al. 1993). Let us study the 
impact of the mechanisms on the spectral lag in the following. 

We first investigate the lag's dependence on the observed peak energy of the vFy spec- 
trum. Ignoring the minute difference caused by various factors, we assume that the following 
relationship holds in any situations concerned in this paper: Ep = l.GYTEop (one can check 
that, for photons emitted from an ejecta moving towards the observer, Ep — 2r£^o,p)- This 
relationship tells us that the observed peak energy is proportional to F or £^o,p when the 
other is fixed. The lag's dependence on the observed peak energy in the case of Ep varying 
with r when fixing Eq^p is similar to those discussed in §2.1. For the sake of comparison, 
we convert the right panel of Fig. 1 to the last panel of Fig. 3 by applying the relationship. 
In addition, we explore the lag's dependence on the observed peak energy in the case of Ep 
varying with £^o,p when fixing F. 
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Figure 2. -Relationships between Lagij,, Lagu, LagQ and the local pulse width for i?o.p=0.75 keV (the left panel), and that 
between Lagn and local pulse width (the right panel) for i?o^p=0.75 keV (the solid line) and 0.15 keV (the dotted line). Here 
we tcike r = 200. Other parameters and symbols axe the same as those adopted in Fig. 1. 

One would find from tlie left panel of Fig. 3 that, similar to Fig. 1, when Ep varies with 
Eq^p for a fixed F, each of the three kinds of lag has a peak lag {Lagp) as well. 

To make clear what causes the relationships shown in Fig. 3, we divide the light curves 
of equation (1) into two parts, with one being contributed from the low-energy portion of 
the rest frame Band function spectrum and the other from the high-energy portion, and 
explore how the peak times, tp, of the light curves are affected by the observed peak energy 
Ep when it shifting through the channel (see Figs. 4 and 5). 

Let us denote the channel concerned by [Ei,E2] (for the first channel of BATSE, it is 
[Ei,E2] = [25,50]keV, while for the third channel, it is [Ei,E2] = [100, 300]A;e\/). From 
Figs. 4 and 5 some conclusions could be reached. For a certain value of the Lorentz factor 
(here we take F = 100), when Ep < Ei, the contributions to the light curves mainly come 
from the rest frame high-energy portion. According to the Doppler effect, this amount of 
energy would mainly come from the area of the fireball surface around the line of sight, i.e., 
9 K, (where 9 is the angle to the line of sight). In other words, energy arising from the 
rest frame high-energy portion (photons with Eq > Eq^p when emitted) emitted from the 
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Figure 3. —Relationships between Lagis, Lagi4, Lagc and Ep for F = 200 (the left panel), and that between Lagn and Ep 
for the two cases of fixing F (the middle panel) and i?o.p (the last panel). In the middle panel, the solid lines from the top to 
the bottom represent F=10, 50, 100, 200, 400, 600, 800, 1000, respectively, and in last panel, the solid lines from the top to 
the bottom stand for £o,p=0.75, 0.15, 0.075 keV, respectively. Other parameters and symbols are the same as those adopted 
in Fig. 1. 



area with ^ would shift into the channel concerned, and thus the peak time tp of the 
corresponding light curves would keep almost unchanged (see Fig. 5). When Ep increases, 
the contributions to the light curves from the low-energy portion become larger since the 
corresponding photons begin to shift into the channel concerned. At the same time, some 
high energy photons emitted from the area corresponding to a larger 9 fall into the channel 
as well, and due to the well-known time delay their contributions to the observed light curves 
cause a larger tp (see Fig. 5). It is the difference of tp of two different channels that leads 
to a lag (see Fig. 1). When Eq^p is fixed, the situation is quite different. In this case, the 
enhancement of Ep comes from the increasing of F. A large Ep corresponds to a large F 
which makes the light curve narrower (i.e., FWHM oc F~^, see Qin et al. 2004; Lu et al. 
2005a ). As a result, tp sharply decreases with Ep (see the upper panels of Fig. 4 and the 
right panel of Fig. 5). Note that, since the influence of the Lorentz factor is very significant 
(as FWHM oc F"^), the difference of tp of two different channels is hardly detectable in the 
right panel of Fig. 5. 

Figs. 6 and 7 show that the lag increases with the low energy spectral index cto of the 




Figure 4. —Plots of the light curves of BATSE channel 3 contributed from the rest frame low-energy portion (the dash line) 
and high-energy portion (the dot line) , and from the whole energy range (the solid line) of a Band function spectrum in cases 
of fixing Bo,p = 0-75keV (the upper panels) and T = 100 (the lower panels). The panels from the left to the right stand for 
Ep=100 keV, 250 keV, 300 keV, respectively. Other parameters are the same as those adopted in Fig. 1. 



rest frame Band spectrum and it decreases with the high energy spectral index /3o of the 
spectrum. The relationships also follow Lagc < Lagi^ < Lagi4 in both cases. Compared with 
that associated with the low energy spectral index, the variation of the lag with respect to 
the increase of /3o is relatively mild, where Lagi^ and Laga are almost unchanged (see Fig. 
7). 

To investigate what causes the relationships shown in Figs. 6 and 7, once more we divide 
the light curves of equation (1) into two parts, as done in Fig. 4, here different ao and /?o 
are adopted (see Figs. 8 and 9). 

Fig. 8 shows that, for BATSE channel 3 and in the case of Ep — 250keV, when Pq is 
fixed, the contribution from the rest frame low-energy portion to the light curve decreases 
slightly with the increasing of ckq, whereas the contribution from the rest frame high-energy 
portion slightly increases with the increasing of ao- The peak time tp of the first channel 
light curve of BATSE rises more rapidly than that of the third channel (see Fig. 8). Note 
that, to shift into the channel observed, photons of the rest frame high-energy portion must 
be those emitted from the area of larger 9, and it must be this that leads to the increasing 
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Figure 5. -Relationships between Ep and the peak time tp of the Ught curves of BATSE channels 1 and 3 in cases of fixing 
r = 200 (the left panel) and -Eo,p = 0-75keV (the right panel). Here we take Rc — D — etc = 0. Other parameters are the same 
as those adopted in Fig. 4. The solid and the dot line stand for BATSE channels 1 and 3, respectively. (Note that the two lines 
in the right panel almost overlap each other.) 

of tp with the increase of Uq. The difference of tp of the two channel hght curves causes the 
lag (see Fig. 6). 

Prom Fig. 9 we find a slightly different situation. For BATSE channel 3 and in the 
case of Ep = 250keV, when ao is fixed, the contribution from the rest frame high-energy 
portion to the light curve is significant only when is close to —2. The peak time tp of the 
first channel light curve is almost unchanged with the increase of /Sq- Only when Ep is low 
enough, the contribution from the rest frame high-energy portion to the light curve could 
become significant. In this situation, the peak time tp of the third channel would rise with 
the increase of (see Fig. 9). This is because that when Ep is low photons of the rest frame 
high-energy portion emitted from the area of larger 6 would shift into the channel observed, 
and then the index would play a role. 

2.3 Spectral lag in the case of a rest frame spectrum varying with time 

It is common that the spectrum parameters of many GRBs are observed to vary with time 
(see Preece et al. 2000). We are curious how the lag's dependence on F and Ate^pwEM 
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Figure 6. -Relationships between Lagis, Lagxi, Lagc and ao for Eo,p=0.75 keV (the left panel), and that between Lagis 
and oq for different -Bo,p (the right panel). Here we take /3o = —2.25 and F = 200. In the right panel, the open circle 

and the open square represent Eo.p = 0.lkeV,0.l5keV , respectively, crosses from the top to the bottom represent -Bo,p = 
0.35keV,0.75keV,lkeV,2keV,3keV, respectively. Other parameters and symbols are the same as those adopted in Fig. 1. 



would be if the rest frame spectrum develops with time. As most of GRBs show "hard-to- 
soft" evolutionary feature (see e.g., Ford et al. 1995; Fishman et al. 1995; Band, David L. 
1997; Preece et al. 1998;), here we consider a simple case where the rest frame spectrum is 
a Band function with its indices and the peak energy Eo^p decreasing with time. Like Qin 
et al. (2005b), we assume a simple evolution of indices ao and /3o and peak energy £^o,p 
which follow ao = -0.5 - k{Te - Te,i)/(re,2 - t^j), /?o -2 - k{T0 - Te,i)/(re,2 - tq^i) and 
loguop = 0.1 - k{r0 - Te^i) / {t0^2 - re,i) (keVhr^) for tq^i < rg < rg^2- For Tg < re,i, ao = -0.5, 

= —2 and logu^p = 0.1 {keVh~^), while for Tg > rg^2, c^o = —0.5 — A;, /5o = — 2 — A; and 
logvQp — 0.1 — k {keVh~^). We take k=0.1, 0.5 respectively (they correspond to different 
rates of decreasing) to investigate the lag's dependence. 

Let us employ local Gaussian pulse (2) to study this issue. We adopt rg^i — %(7 -\- Te^min 

and Te,2 ^11(7 + Tg^min and Te,o = lOcT -\- Tg^rnin and Tg^rnin = 0. 

Figs. 10 and 11 show the dependence of lag on the Lorentz factor and the local pulse 
width in the case when parameters of the rest frame spectrum vary with time. We find that, 
when the spectral parameters decrease with time, the lags are larger than that associated 
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Figure 7. -Relationships between Logia, Lagu, Laga and /3o for Eo,p=0.75 keV (the left panel), and that between Lagis 
and Po for different £^o,p (the right panel). Here we taJje ao = —1 and T = 200. Other parameters and symbols are the same 
as those adopted in Fig. 6. 



with the case when the spectral parameters are fixed. The larger the decreasing speed (say, 
k=0.5), the more obvious the effect. 

Kocevski et al. (2003a) argued that the observed lag is the direct result of spectral 
evolution and found the observed lag increases with the peak energy's decay rate. Motivating 
by this, we take ao = —0.5, /5o = —2, and allow E^p varying with time. Illustrated in Fig. 12 
is the relationship between the Lagis and the rest frame peak energy decaying rate, k, which 
is in agreement with the observed result found by Kocevski et al. (2003a). The dependence 
of lag on the varying speed k in the case of the rest frame spectrum parameters ao and (3o 
varying wit time is also shown in Fig. 12. We find it common that the lag increases with the 
increasing of k. The increasing of lag in the case when the rest frame peak energy decreasing 
with time is much more obvious than that in the case when indexes ao and (3o decreasing 
with time. 



2.4 Lag's dependence on energy 

In a previous study, Wu & Fenimore (2000) detected no larger lags in lower energy bands. 
However, we find from the above analysis that the following relation holds in almost all the 




Figure 8. -Plots of the light curves of BATSE channel 3 contributed from different portions of the rest frame Band function 
spectrum for Ep = 250keV (the first three panels) and the relationship between the peak time of the light curves of BATSE 
channel 3 (the solid lines) and 1 (the dot lines) and ao for different Ep (the last panel). Here we take Po = —2.25 and F = 200. 
In the first three panels, the dot, dash and solid lines represent the light curves contributed from the high-energy portion, the 
low-energy portion and the whole energy range of the adopted spectra, respectively. The solid lines (or the dot lines) from the 
bottom to the top in the last panel stand for Ep = 50, 100, 250, 300, 400 keV, respectively. The Other parameters are the same 
as those adopted in Fig. 5. 

cases concerned: Lagc < Lagi^ < Lagu. Motivated by this, we investigate in the following 
in much detailed the time lag's dependence on energy. We adopt a typical Band function 
spectrum with — —1, /Sq = —2.25 and Eq^^ — Q.lbkeV as the rest frame radiation form 
to study this issue. 

There are three instruments in the Swift telescope pay load, and two of them could detect 
high energy photons with the XRT ranging from 0.2 keV to 10 keV and the BAT ranging from 
15 keV to 150 keV (http:/ / heasarc.g s f c.nasa.gov / docs / swi ft / aboutsWi ft /) . We consider a 
wide band ranging from 0.2 keV to 8000 keV which covers both the Swift and BATSE bands. 
This band is divided into the following uniformly ranging channels (where E2 — 2Ei): [Ei, 
E2MO.2, 0.4] keV, [0.5, 1] keV, [1, 2] keV, [2, 4] keV, [5, 10] keV, [10, 20] keV, [20, 40] keV, 
[50, 100] keV, [100, 200] keV, [200, 400] keV, [500, 1000] keV, [1000, 2000] keV, [2000, 4000] 
keV and [4000, 8000] keV. 

We measure the lags between the first channel and any of the other channels. The 
relationship between the lag and the corresponding energy E (here E denotes the lower 
energy limit of the corresponding high-energy channel) is presented in Fig. 13. It shows 




Figure 9. -Plots of the light curves of BATSE channel 3 contributed from different portions of the rest frame Band function 
spectrum for Ep = 250keV (the first three panels) and the relationship between the peak time of the light curves of the light 
curves of BATSE channel 3 (the solid lines) and 1 (the dot lines) and /3o for different Ep (the last panel). Here we take ao = — 1- 
Other parameters and symbols are the same as those adopted in Fig. 8. 



that, for any given value of T, the lags increase with energy following the law of lag oc E 
and then saturate at a certain energy (we call Es, where subscript "s" represents the word 
"saturate"). This might explain why Wu & Fenimore (2000) found that GRBs do not show 
larger lag at lower energy. We notice that Eg takes place behind the corresponding peak 
energy Ep (where the relation Ep — 1.67r£'o,p is adopted). For different F, Es increases with 
Ep following the law of Eg oc Ep (see the right panel of Fig. 13), and the relation Eg oc Ep 
holds for all -Eo,p- 

As shown in §2.2, when E > Ep, the contributions to the corresponding light curve 
largely come from the high-energy portion of the rest frame spectrum, which makes the 
peak time of the light curve less change, and thus the lag would saturate. The contrary is 
true. When E <^ Ep, the contributions come only from the low-energy portion of the rest 
frame spectrum, emitted from the whole fireball surface, i.e., < ^ < tt/2. Thus, there 
would be almost no lag between the light curves of such two energy channels (see the two 
filled circles in the bottom of the left panel of Fig 13). 
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Figure 10. -Relationship between Lagn and F associated with the situation when the rest frame spectrum varying with time 
in the case of k = 0.1 (the left panel) and k = 0.5 (the right panel). For the sake of comparison, the relationship deduced 
from the constant rest frame spectrum with the corresponding upper and lower limits of the indexes are also presented. The 
implications of the lines are as follows: v-cons (ao and /3o varying with time, loguQp = 0.1 keVh~^); ab-cons (ao = —0.5 and 
(3o = —2.0, uop varying with time); abv {ao, 0o and i^op varying with time); 52001 (ao = —0.5, po = —2.0 and logvQp = 0.1 
keVh-^); 62100 (ao = -0.6, /3o = -2.1 and loguop = 0.0 keVh-'^); 12504 (ao = -1.0, /3o = -2.5 and logvop = -0.4 keVh-^). 
Here we take At0^p]^f{M = lO^s. Other parameters are the same as those adopted in Fig. 1. 

2.5 Lag's dependence on the opening angle of uniform jets 

According to the relativistic fireball model, emissions from a spherically expanding shell and 
a jet would be rather similar to each other as long as we are along the jet's axis and the 
Lorentz factor F of the relativistic motion satisfies F"^ < 9j, because the matter does not 
have enough time (in its own rest frame) to expand sideways at such a situation (Piran 1995, 
1999). Let us study the dependence of lag on the radiated area confined by the opening angle 
of Omax — p/^-i where p is a parameter describing the width of the opening angle of uniform 
jets. Shown in Fig. 14 is the relationship between Lagi^ and p. 

We find from Fig. 14 that Lagiz increases with the opening angle within Omax < O.GF"^ 
following the law of Lagi^ oc p^-^, and beyond this range, it slightly increases with the angle, 
and finally saturates at p = 0.8 for any F. The saturate lag (denoted by Lags) decreases 
with the increasing of F following the law of Lags oc F~^'^ (see the right panel of Fig. 14). 
This is an outcome of the beaming effect. It suggests that there is a limited value of the lag 




Figure 11. —Relationship between Lagyi, and the local pulse width Atg fwhm associated with the situation when the rest 
frame spectrum varying with time in the case of = 0.1 (the left panel) and k = 0.5 (the right panel). For the sake of 
comparison, the relationship deduced from the constant rest frame spectrum with the corresponding upper and lower limits 
of the indexes are also presented. Implications of the lines are the same as those in Fig. 10. Here we taJje F = 100. Other 
parameters are the same as those adopted in Fig. 1. 

when different opening angles of uniform jets are concerned, which is the same as that in 
the case of spherical fireballs. 

2.6 Lag's dependence on the radius of fireball 

Shen et al. (2005) found that the lag is independent of the radius of fireballs based on their 
model (see Fig. 8 in their paper). We perform the analysis with their formula when raising 
a higher precision of calculation (the result is the same when we adopt equation [1]). Here 
we fix the local pulse width by adopting Atg = 10s. The result is presented in Fig. 15. 

Fig. 15 shows that Lagi^ decreases with the radius of fireballs following the law of 
Lagis when the local pulse width is fixed. The lag depends on the radius of fireballs. 

Comparing our result with that obtained in Shen et al. (2005) we find that the conclusions 
are entirely different. We suspect that three factors might be responsible to this difference. 
The first is the precision of calculation. The second is the range concerned. Our calculation 
spans about three magnitudes of the fireball radius while Shen et al. (2005) considered only 
two magnitudes. The third is that the dependence of Lagis on the radius is weak: the lag 
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Figure 12. -Relationship between Lagis and the varying speed k of the parameters of the rest frame spectrum. The implications 
of the lines are as follows: v-cons («o and 00 varying with time, logEop = 0.1 keVh~^); ab-cons (ao = —0.5 and /3o = —2, Eop 
varying with time); abv (ao, ft) and uyp varying with time). Here we taJie F = 200 and Atg^p-^YHM = lO^s. Other parameters 
are the same as those adopted in Fig. 1. 



would be only about two times smaller when the radius becomes three magnitudes larger. 
We believe that the third is the key factor accounting for the difference. 

The conclusion of Shen et al. (2005) that the lag is independent of the radius of fireballs 
is a puzzle since a large radius seems to cause a larger distance between the area with ^ = 
and that with > on the fireball surface and this should lead to a larger lag. Our conclusion 
that the lag decrease slightly with the increasing of the radius makes the puzzle worse. The 
mechanism accounting for this is currently unclear. We suspect that this dependence might 
rely on what kind of photons dominate the peaks of the lower and higher energy channels 
respectively and where they are emitted from. In creating the dependence, other parameters 
such as the rest frame peak energy and the Lorentz factor might be at work. Therefore, this 
might not be answerable if only a simple mechanism is concerned (to find a answer to this, 
a detailed investigation should be made, which is beyond the scope of this paper) . 
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Figure 13. —Relationship between the lag and the lower energy limit E of the corresponding high energy channel for different 
r (the left panel) and that between Es and Ep (the right panel), where the relationship Ep = ISJVEo.p is adopted and we 
take oq = ~1) A = —2.25 and i?o,p = 0.75keV. Solid lines from the top to the bottom in the left panel stand for r=10, 50, 
100, 200, 500, 1000, 2000, respectively. The solid line in the right panel is the linear fit one of the data. Other parameters are 
the same as those adopted in Fig. 1. 

3 DIFFERENT CHARACTERISTICS SHOWN IN THE TWO CASES OF 
FIXING THE INTRINSIC PULSE WIDTH AND THE LOCAL PULSE 
WIDTH 

Shen et al. (2005) explored the lag's dependence on the Lorentz factor and found Lag oc F"^ 
when fixing the intrinsic pulse width (the co-moving pulse width) and the fireball radius. 
However, in this paper we find Lag oc F^^ instead, when fixing the local pulse width and the 
fireball radius based on both models of Shen et al. (2005) and Qin et al. (2004). According 
to the theory of special relativity, the two results are identical, as the local timescale is F 
times of the corresponding intrinsic timescale. Some interesting questions arise accordingly. 
When the fireball radius is fixed (as considered in Shen et al. 2005 and in this paper), would 
the two approaches always lead to identical results? Could they reveal the same property of 
the curvature effect? In which case should which one be preferred? 

As pointed out and illustrated in Qin et al. (2004), the profiles of light curves are not 
affected by the Lorentz factor when fixing the relative local pulse width {cAtg/Rc, or Arg as 




Figure 14. -Relationship between Lagia and p for different F (the left panel), and that between lags and F (the right panel), 
where we adopt a typical Band function spectrum with ao = —1, /3o = —2.25 and Bo,p = 0.75keV as the rest frame radiation 
form. In left panel, solid lines from the top to the bottom stand for r=100, 300, 500, 700, 1000, respectively. Other parameters 
are the same as those adopted in Fig. 1. 

used in Qin et al. 2004). When the fireball radius is fixed, a fixed relative local pulse width 
would correspond to a fixed local pulse width (that is, a certain Arg leads to a certain Atg in 
this situation, where Rc is fixed). Thus, in this situation, fixing the local pulse width would 
not lead to difi^erent profiles of light curves when different Lorentz factors are considered. 
However, the feature of the profile will change when fixing only the intrinsic pulse width. 
This is because that for a ceratin intrinsic pulse width, the corresponding local pulse width 
would take different values when different Lorentz factors are considered (recall that the 
local timescale is F times of the corresponding intrinsic timescale), and this would lead to 
different relative local pulse width since Rc is fixed. The corresponding profiles of light curves 
would change since they are sensitive to the relative local pulse width (see Qin et al. 2004). 

Fig. 16 illustrates the different infiuences of the Lorentz factor on the profiles of the light 
curves calculated with equation (5) in Shen et al. (2005) in the two cases of fixing the local 
pulse width and the intrinsic pulse width. 

Kocevski et al. (2003b) and Ryde et al. (2003) found that there is a linear relationship 
between the rise width, FWHMj. and the full width, FWHM, of gamma-ray pulses. Lu et 
al. (2005a) explained the relationship based on the formula presented in Qin et al. (2004), 




Figure 15. —Relationship between Lagis and the radius of the fireball. Here we take T = 100, Eo,p = 0.75keV and 
^tg^p\^HM = lOs- The solid line is a linear fit to the data. Other parameters are the same as those adopted in Fig. 1. 



and found that there exists a dead hne in the FWHMr — FWHM plane when taking the 
relative local pulse width Atq^fwrm > 1- Here we repeat the same work as illustrated in 
Fig. 2 of Lu et al. (2005a) based on the Shen et al. (2005) formula. As shown in Fig. 17, 
we find that there indeed exists the dead line predicted before in the FWHMr — FWHM 
plane when fixing the local pulse width (here we take td > lOOOs and = 3 x lO^^cm 
which corresponds to Atq^fwhm > !)• However, when fixing the intrinsic pulse width, one 
cannot find a dead line in the plane (see Fig. 17). This indicates that the features in the 
FWHMr — FWHM plane are different in the case of fixing the local pulse width and in 
that fixing the intrinsic pulse width. One could obtain the local pulse width td by measuring 
the location of the corresponding data (to observe in what lines they belong to) from the 
left panel of Fig. 17 as long as Atq^fwhm < 1- But one could not measure the intrinsic pulse 
width from the right panel of the figure. This indicates that the former is superior to the 
latter when one attempts to make use of the FWHMr — FWHM plane. 




Figure 16. -Plots of the light curves of BATSE channel 3 resulting from symmetric Gaussian pulses when fixing the local 
pulse width td = 100s (the left panel) and the intrinsic pulse width td' = 10s (the right panel) for different F. Here we take 
a = —1, /3 = —2.25, Rc = 5 x 10^'*cm and Ep = 350keV (where the relationship Ep = LGlFEo p is adopted). The solid, dash 
and dot lines stand for F = 10, 100, 1000 respectively (the three lines in the left panel overlap each other). Note that, for the 
sake of comparison, we normalize the light curves and re-scale the variable so that the peak count rate is located at T' = and 
the FWHM of the decay portion is located at T' = 0.2. 



4 DISCUSSION AND CONCLUSIONS 

Discussed in Shen et al. (2005) and the above analysis, only intrinsic or local Gaussian 
pulses are involved. One might ask if the conclusions hold only in the selected pulse form. 
Thus, we replace the local Gaussian pulse with other forms of local pulses when applying 
equation (1) and the corresponding formula of Shen et al. (2005). Local pulses (55), (81), 
(82), (83), (85), (86) and (87) presented in Qin et al. (2004) are studied. We find that the 
lags resulting from a local pulse without a decaying phase are too small to be noticed. For 
example, laglSt < 10~^s when F > 100 and Rc — 3 x lO^^cm. This is in well agreement 
with what discovered by Shen et al. (2005). Now it is clear that, in the mechanism of the 
curvature effect, it is the decaying phase of a local pulse that contributes to the observed 
lags. 

A recent study revealed that (see Lu et al. 2005b), the rise phase of light curves will 
always be dominated by the emission from 6' ~ and within the rise phase of the local 
pulse. It is illustrated in Fig. 18, where a local pulse without a decaying phase is considered. 
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Figure 17. Relationships between FWHMr and F\V HAl for tiic light curves of BATSE channel 3 resulting from symmetric 
Gaussian pulses associated with the case of fixing the local pulse width (the left panel) and the case of fixing the intrinsic 
pulse width (the right panel). In the left panel, solid lines from the bottom to the top represent td = 10s, 100s, 1000s, 10000s, 
respectively (note that the last two solid lines overlap each other), and in the right panel, the corresponding lines stand for 
td' = 0.1s, Is, 10s, 100s, respectively. Here we take Rc = 3 x 10^^ cm and other parameters are the same as those adopted in 
Fig. 16. 

This indicates that there would be no peak lags for the light curves of different energy 
channels associated with those local pulses without a decaying phase. 

Except the Band function spectrum, the Comptonized spectrum and the broken power 
law spectrum were employed to fit the spectral data in Preece et al. (2000). To make sure if 
the dependence of lags discussed above would be altered when other rest frame spectra are 
considered, we employ the Comptonized spectral model as the rest frame radiation form to 
investigate the lag's dependence concerned above. The Comptonized spectrum is written as 
Ii, oc i^o~fi°'°''^ exp{— , where q;o,c and z/q.c are constants. Here, typical values, «o,c = —0.6 
(Schaefer et al. 1994) and uo,c = O-bbkeVh'^ (see Qin 2002 Table 2) are adopted. Our 
calculation shows that the characteristics of the lag's dependence are the same as what 
revealed above, when the rest frame Band function spectrum is replaced with the two other 
forms. The only difference is the magnitude of the lags. For example, when the rest frame 
Comptonized spectral model is adopted, we obtain a slightly larger value of Lagi^ (about 
70% larger) when we redraw Fig. 2. 

Of the sample of 41 source presented in Friedman et al. (2005), we find that the range of 
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Figure 18. Plots of the normalized light curves of BATSE channel 3 arising from a local linear rise pulse with Arg^pwHM = 1 
(a local pulse without a decaying phase) emitted from different areas of the fireball surface. Here we take T = 200 and adopt a 

Band function rest frame radiation form with cto = —1, /3o = —2.25 and UQp = 0.75keVh^^ . The solid lines from the right to 
the left stand for Smaa;=T/2, l/F, 0.5/r and O.l/F, respectively (note that the rising phase as well as the peak of the counts 
rates for different light curves overlap each other). Other parameters are the same as those adopted in Fig. 1. 

6j is {0.03107, 0.57072 } rad, and its average value is 0.15925 rad. It shows that 6j ^ F"^ if 
r > 100. As discussed in §2.5, one could find that for GRBs coming from a uniform jet with 
such opening angles, if F > 100 is adopted, the difference of lags from that arising from a 
spherical surface would not be distinguishable. 

Several authors showed that the CCF lags between BATSE channels 1 and 3 tend to 
concentrate near < 200ms (see Norris et al. 1996, 2000; Chen et al. 2005). A recent mea- 
surement carried by Norris et al. (2005) for BATSE wide pulses showed that the long lag 
range is 1.0 < lag < 4.2 s. It is generally believed that the gamma-rays arise from internal 
shocks at a distance of ~ 10^^ — 10^^ cm from the initial source (Piran 2005). Within 
the generally accepted ranges of Rc — 10^^ — 10^^ cm and F > 100, the observed lags could 
be accounted for by the curvature effect when different values of parameters are adopted. 
For example, one gets Lagi^ = 300 ms from Fig. 1 when taking F = 100; Lagi^ = 100ms 
from Fig. 7 when taking ao = —1, /3o = —2.25 and £'o,p = 0.75keV. According to Fig. 
2, Lagis — 4.6 x 10~'^ Ate^pwRM s. One could obtain a large lag, i.e., Lagis — 4.6 s when 
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taking Ato^pwHM = 10''s. In this situation, the FWHM of the corresponding hght curve 
determined by equation (1) is about 75 s when taking F — 500, which is very close to the 
observed data of the wide-pulse bursts studied in Norris et al. (2005), i.,e., 1.0 < lag < 4.2 
s and 6s < the width of channel 1 < 100s (see Fig. 6 in that paper). All these show that 
merely the curvature effect could produce the observed lags. A larger lag requires a wider 
local pulse. 

We thus come to the following conclusions: a) lag oc F^*^ with e > 2 within different 
energy channels, and the relationship would different for different rest frame peak energy 
Eop] b) lag is proportional to the local pulse width and the FWHM of the observed light 
curves, and the local pulse width plays a more important role in producing a large lag when 
compared with other lag's dependent parameters; c) a large lag requires a large ccq and a 
small /3o as well as a large Eop] d) when the rest frame spectrum varies with time, the lag 
would become larger; e) lag decreases with the increasing of Re, f) lag oc E within a certain 
energy range for a given F; g) lag is proportional to the opening angle of uniform jets when 

e^ax < 0.6F-1. 
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